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Abstract 

Based on the Tsalhs entropy, the nonextensive 
thermodynamic properties are studied as a q- 
deformation of classical statistical results using 
only probabilistic methods and straightforward cal- 
culations. It is shown that the constant in the Tsal- 
lis entropy depends on the deformation parameter 
and must be redefined to recover the usual ther- 
modynamic relations. The notions of variance and 
covariance are generalised. A partial derivative for- 
mula of the entropy is established. It verifies impor- 
tant relations from which most of the nonextensive 
thermodynamics relations can be recovered. This 
leads to a new proof of a highly nontrivial con- 
clusion that thermodynamic relations and Maxwell 
relations in non extensive thermodynamic have the 
same forms as those in ordinary nonextensive sta- 
tistical mechanics. Theoretical results are applied 
to ideal gas. The case of fermions and bosons sys- 
tems with fractal distributions is also considered. 

PACS : 05.20.-y, 05.70.-a, 05.30.-d. 



tems involve long-range interactions, statistical me- 
chanics presents serious difficulties due to non- 
markovianity of stochastic processes. Non-extensive 
thermodynamics may correctly cover some of the 
known anomalies. The canonical ensemble of non- 
extensive thermodynamics for quantum mechanical 
systems with finite number of degrees of freedom 
has been studied in Ref . |^ . It has been shown that 
most of thermodynamic relations are preserved |^ , 
if introducing the entropic index q, we replace the 
usual means by g-expectations. However, invariant 
relations involve Bolzman-Gibbs constant in the 
canonical ensemble and are related to Tsallis con- 
stant in the grand canonical ensemble. The canon- 
ical distribution in non-extensive thermodynamic 
is given by optimal Lagrange multipliers and is dis- 
cussed in Ref. Q . Parameter differentiation of frac- 
tional powers of operators has been developed in 
Ref. 1^, ^ . Nonextensive statistical mechanics and 
thermodynamics, as a new stream in the research 
of the foundations of statistical mechanics, is pre- 
sented in a recent publication ||ic|l . 



1 Introduction 

The purpose of this work is to discuss the ther- 
modynamic relations problem in Tsallis thermo- 
statistics recently presented in Ref. [|l|. It has 
been pointed out that (/-thermostatics could be 
understood in terms of g-deformations. The Tsal- 
lis' entropy proposed by Constantino Tsallis ten 
years ago |Q has found many applications in var- 
ious domains of physics [Q, as the turbulence in 
electron plasmas, phonon-electron thermalization 
in ion bombarded solids and flux of solar neu- 
trinos. There are also applications to the mech- 
anisms of anomalous diffusion B. When sys- 



2 Tsallis Entropy 

The Tsallis' entropy is defined by the formula 

Sg ^ k J plngip-^) dX (1) 

where p is the distribution of probability, A is the 
Lebesgue measure on the phase space and fc is a 
constant, called the Tsallis constant. In this ex- 
pression the usual logarithmic function is replaced 
by the q-logarithmic function 



which is the inverse function of the g-exponential 
function 



e,{x) = (1 + (1 - 



(3) 



The problem posed by unnormahzed expectation is 
studied in Ref. Introducing the q-Gibbs mea- 
sure (Ref. [Q) in the canonical ensemble 



defined if (1 + (1 — g)a:) > and which is prolonged 
by continuity with eq{x) — otherwise. When q 
1, the ^-exponential and the g-logarithmic func- 
tions converge to the ordinary exponential and log- 
arithmic functions. The entropy of Tsallis reduces 
to the Bolzmann-Gibbs formula in the extensive 
limit. The pseudo-additivity of the g-logarithmic 
function 

lnq{xy) = ln,(a::) + lnq{y) + {1 - q) \ng{x) lng{y) 

(4) 

leads to the nonextensivity of the Tsallis entropy. 
For two independent subsystems A and B of a total 
system, in thermal contact with each other, it has 
been shown that 

Sq{pA Pb) = SqipA) + Sq{pB) + (1 ^ q)kSq{pA)Sq 

(5) 

The following property of the q-logarithmic func- 
tion 



P = 



Z„ 



(12) 



where H is the Hamitonian of the system and /? 
the Lagrange multiplier (l/kT). The expression of 
the density of probability has been discussed in Ref. 
[p^.It is shown in Ref. [|l5| that, for a canonical dis- 
tribution of a system in contact with heat bath, the 
entropic index q is related to the particle number 
N and characterizes the bath. The Tsallis constant 
k is choosen as 



g 

CqksT 



(13) 



where ks is the Boltzamm-Gibbs constant. The q- 
■^jf^^ition function Zq is given by the normalisation 
of the probability density 



1 



Zq = J eq{-l3H) dX 



(14) 



ln,(-) = -x^-Hnq{x) 



(6) 



gives the usual definition of the Tsallis entropy (TT| 

Sq = -k I pnUqip) dX (7) 

which is also defined by the equivalent formula 



S„ = 



1 - p" dX 



(8) 



If we define the unnormalized g-expectation of an 
observable A by 



\A]q^J Ap^ dX 



and let Cq be the constant 



p« dX 



(9) 



(10) 



The measure fi defined relatively to the Lebesgue 
measure by dp, = (p"^ /cq) dX is a measure of proba- 
bility. In the extensive limit g — > 1, the q-Gibbs 
measure converges to the ordinary Boltzmann- 
Gibbs density. The substitution of the g-additivity 
Eq. (^) in the Tsallis entropy Eq. (|l|) gives 

Sq = -kJ {-l3H + \nq{Z-^) - (1 - q)l3H\nq{Z-^))p'' dX 

(15) 

and using Eq. (||) 

Sq^k j{f3H + {l-{l- q)l3H)Zl-^ \^<i{Zq))p'i dX 

(16) 

we find that the Tsallis entropy 



the g-expectation (normalised) could be defined by 
= = 1 / Ap^ dA (11) 



Sq = k j {fiHp'i + \nq{Zq)p) dX (17) 

leads to the following formula 

Sq = kP \H^^ + k\r,q{Zq) 
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Introducing the q-internal energy 
and the g-free energy 

-1 



(18) 



F, = —ln,iZ,) (19) 



one finds 



Fq = CqJJq - TSq 



(20) 



This relation is the q-generahsation of the ordinary 
equation of 5 = U/T + kln{Z). 



3 Mathematical Formalism 

In this section, we introduce five important rela- 
tions for the nonextensive mathematical formalism. 
The generalised internal energy Uq defined as the 
^-expectation of the Hamitonian is given by 



1 dlnqjZq) 
Cq d/3 



(21) 



The proof of this equation is easy to compute 

d\nq{Zq) _ d\nq{Zq)dZq _-ldZ 



22) 



5/3 dZq df3 Z| a, 

- ^|/e.Mi^)^A (23) 

Using the following expression of the derivative of 
the (/-exponential function with respect to (3 



deq{-l3H) 
dp 



-Heli-pH) (24) 



leads to 



dp 



^q 



J Hp" d\=\H}^^CqU{,2Q) 



Defining the g-movariance of two operators A and 
Bby 

Movq{A,B) = \AB'\^^_, - \A]^ \B]^ (27) 



the equation (Eq. gives the g-movariance of A 
and H in function of the derivative of the g-mean 
of A 



dJAl 
dp 



OA 
'dp 



-qZ^-Hlovq{A,H) (28) 



The movariance converges to the usual covariance 
in the extensive limit g — > 1. Proof of the equation 
(Eq. H) is 



d\A\^ 
dp 



dA 
'dp 

dA 

'dp 



p'i dX 



Ap- 



^djf_dp_ 
dp dp 
^_,d (eq{-PH) 



(29) 



dp 



Za 



The derivative of the probability density with re- 
spect to the Lagrange multiplier is computed as 



5p 

dp 



1 deq{-pH) eqi^pH) dZq 



Za 



dp 



^q 



dp 



-H 



eq{-PH) 



Substituing the previous equation in Eq. 
using that Zq is a constant leads to 



(31) 

-PHim 

J) and 



d_\Al 

dp 



dA 
'dp 



(33) 



qz^-^ j AHp^"-^ dX + - 

iZ^^ J Ap'i (^J Hp'' dX^ d^4) 
qZ"^-\\AH^,^_,~\A^^\HM^5) 



-qZl-Hlovq{A,H) 



(36) 



If the operator A does not depend explicitly of /3, 
we find that the g-movariance is proportionnal to 
the derivative of the g-expectation of A with respect 
to p. We define the g-mariance by 



Marq{A) = \A^] 



2q~l 



(37) 



This function converges to the usual variance when 
q tends to 1. Using the equation of the generalised 
energy Eq. (^), we see that the q-mariance of the 
hamiltonian is given by 



Marq{H) 



1 dHnq{Zq) 



qzr' 



dp^ 



(38) 
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Applying Cauchy-Schwarz inequality 



The derivative of the partition function with re- 
spect to the parameter is given by 



fg dx 



< 



f dX 



9^ dX 



(39) 

^/■^ proves that 



da 



-a 



— e^(-/3i7) dX 



with / Ap'J-i/^^-g-i and g = 
Marq{A) is always positive. In the extensive limit Let A^^^^ be the quantity 
when q tends to 1, the variance of the hamitonian 
Var{H) = 9^ \n{Z)/dl3'^ is positive. Thus the inter- 
nal energy as a function of the Lagrange multiplier 
/3 is a convex function. 



Ag,Qv4 — 



da 



OA 
da 



(47) 



(48) 



Substituing Eq. (^) and (|47|) in the previous equa- 
We suppose now that the Hamiltonian H and the tion and using that Zq is a constant, we get 



probability density function p depend on a param- 
eter a distinct of (3. 



d\nq{Zq) 
da 



= /3 



dH 
da 



,9-1. 



d feqi~l3H) 



da 



Z„ 



dX 



(40) 



-qZ^-^P j A^-^p^^-^ dX + -- 



Remark that this equation is also valid for a equal P 
(by Eq. (pi])). Proof of this equation is easy to 
compute. 



dlnqjZq 
da 



-IdZq 
Zq da 



-qZ^'P 



eq{-(3H) dX 
(41) 



A 



Ap" 
dH 



dH 



da 



2q-l 



p" dX 



(49) 
(50) 
dX (51) 



dH 
da 



dH 

-qZl-HloVq{A,—) 



(53) 



Using the expression of the derivative of the q- 
exponential function with respect to a 



deq{-l3H) 



leads to 

d\nq{Zq) 



dp 



da 



The following equation characterizes the derivation 
of the Tsallis entropy with respect to the a param- 
eter 



dH 

da 



p'i dX^p 



dH 
da 



(42) 



(43) 



da 



da 



dH 
da 



dH^ 



= qkp^Zl-HloVq{H,—) 



(54) 
(55) 



The next fundamental equation gives the expres- the extensive limit when q tends to 1, this equa- 
sion of the movariance of an observable and the ^i^^ ^^^ws that the derivation of the usual entropy 



derivation of the hamiltonian with respect to a 
given parameter 



is proportional to the variance of the hamiltonian 



da 



dA 
da 



dH^ 



dS 1 1 (d{H) /dH\ 

— = -Var{H) = - —— -( — ) 



- -qpZ^-'MoVqiA, — ) (44) 



da T 



T 



dV \dV / 



(56) 



The proof of this equation is 



da 



dA 
da 
dA 



p" dX 



da 



dpi dp 

A-^-^ dX 

dp da 



If the parameter a is the volume V and the pressure 
defined by F = — {dH/dV) , this equation is the 
usual thermodynamic relation TdS — dU + PdV . 



we define the g-pressure by 



dp 



d (eq{-PH) 



(46^ 



P =-/^\ 1 d\nq{Zq) 



\dV I 



PCq 



dV 



(57) 
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and write the derivative of the TsaUis's entropy 
with respect to the volume 



^ - -kBc 

dv " 



dUa 



leads to the ordinary equation 

dUq = -CqTdSq - PqdV 

Define the following q-quantities 

dQq^ djH)^ /dH\ 
dV dV \dV / ^ 

^^_P ^/^\ 
dV " \dV / 



(59) 



(60) 



(61) 



(62) 



leads to the usual equation in the q-world 

dUq ^ dQq + dWq = dQq - PqdV (63) 

Derivating Eq. (0) with respect to parameter a 



dSq k dcq 



da 1 ~ q da 
Substituing the following equation 

1-9, 



Cq = 1 



k ^« 



(64) 



(65) 



and taking A = 1 in Eq. (Q) leads to the differen- 
tial equation Allows by the Tsallis entropy. 



da da " 

9 .or.a-l(^\H^,^^l 



q 



-kpz-i- 



da 



da 



(66) 
« '''67) 



4 Ideal Gas 



The classical ideal gas in D-dimensional space has 
been exhaustively discussed in Ref. [|l6|. Correla- 
tions were analyzed in Ref. . We summarize 
some results in the case < g < 1 and show by 
straightforward calculation how the partition func- 
tion is computed. The hamiltonian is 



N 



^ 2m 



(68) 



where m is the particle mass, A'' the particle num- 
ber and Pi the momentum of the ith particle. The 
partition function is given by 



{l-a{pl + ...+pI))^-^ 



(69) 



where a is the positive constant a = (1 — q)j3/2m 
and SI is domain of integration such that 1 — a{p\ + 

...+pI)>q. 



z„ = 



V 



N 



1 J^ll^DN 



il-aipl + ...+p%))^ d°p. 



(70) 



and 



(71) 

Let X be ar^ . Using the fact that the beta distribu- 
tion is a density of probability for two parameters 
u and V greatest or equal 1, that is 

After a simple calculation, we find the following 
expression of the (/-partition function 



f 2nrn \ ^^^^ r(f5f) 

m [{i-q)i3h^) r(f5f + ^) 



(73) 



The same calculation for the constant Cq leads to 



Cq = 



l-q ^ '^l-g ^ 2 ) 



(74) 



^ V l-g 2 I ^ y l-q' 

Using the following property of the gamma function 

r(z + i) = zr(z) (75) 

we find 

(76) 



Cq^Z^-^il + il-q)^) 



The generalised energy given by Eq. ( pi|) is com- 
puted in the same way 



'~ 2/3 Cq 



(77) 
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Using the expression of the Lagrange mutipher 
given by Eq. (p9|), we obtain the foUowing expres- 
sion of the internal energy which is independant of 
the entropic index. 



(78) 



The Tsalhs and Boltzmann-Gibbs constants are re- 
lated to the equation 

k^CgZ^-^kB= (l-il-q)^)kB (79) 



(80) 



The equation of state 

1 dZ, _N Z\-'i 
" PcgZl dV V PCq 

leads to the usual expression 



PqV = NkeT 



(81) 



In the extensive limit when q tends to 1, we recover 
the equation of ideal gas. Remark that if we replace 
the hamiltonian by the same expression minus the 
generalized energy, the results are unchanged | p^ . 
The g-mariance of the hamitonian 



Marq(H) = -keTUn 

q 



(82) 



converges in the extensive limit to the variance of 
the hamiltonian DNk^T^ /2. Simple calculations 
give the movariance of the g-pressure 

MoVq{H, Pq) = l—lB^P^CqZl-" (83) 



and 



q 2/3 



^2-29 (84) 



and the ratio 



Mov,{H,P,) 
Mov.iH,^) 



= (9-1) 



k DN 



(85) 



5 Quantum Free Particle 

In this section, a free particle confined in a 3- 
dimensional box is considered. Denoting ii, L2, 



L3 the sides of the box, it is well known that the 
energy eigenvalues are given by 



E = 



f n\ "3 



8m \Ll Ll Ll 



(86) 



where the quantum numbers are non-negative inte- 
gers rii = 0,1,2, ...The partition function is given 
by 

Z,= e,(-/3i?)= (l-(l-'?)i^)^ 

ni,»2,n3 "1,12, "3 

(87) 

Let \t be the translation temperature 



SmkL"^ 



where L Li is the length of a side. For sufficiently 
large T ^ Xt, the sum over rii can be replaced by 
an integral 



Z„ 



1 



(1 - (1 - q)PE) 1-" Inint)dnidn2dn3 

(89) 



where r2 = {1 - (1 - q)f3E > 0} and 1a{x) = 1 if 
X G A and otherwise. For a volume V = L1L2L3 
and a parameter b — /3(1 — q)h'^/8m the change of 
variables 

Y o„3/2 fl/Vb 

Z„ = l^±^ i {l-br^)^.r^dr (90) 



' 8 r(3/2) Jo 
could be written as a beta integral 

Z„ = — ,„,„„,. / {l-x) — x^dr (91) 



63/2r(3/2) 7o 



thus 



27r 



r(f5f) 



this is the classical formula with = 1 and D = ?>. 

6 Fermions and Bosons Distri- 
butions 

Thermodynamics on quantum groups and thermo- 
dynamics of boson and fermion systems with fractal 



6 



distributions have been studied recently [|T9[ ^ . 
In the grand canonical ensemble, the grand parti- 
tion function is the sum under the states j 

e,=j2{i-p{i-qm-fi))^ (93) 

where fi is the chemical potential. The density of 
states in the phase space is given by 

In the thermodynamic limit, the sum is replaced by 
a continous integration, thus the average is defined 

by 

(^>, =E ^^Pj ^ / AD{E)r{E)dE (95) 
■Jo 



where / is the function 
fiE) 



-In (96) 

{1 - (3{1 ^ q){E - fi)) — 

The change of variable E — / (3{q — 1) gives 



(98) 



(99) 



- JO 
The average total number of particles 

1 



, {l-P{l-q){E,~^i))- 
in the thermodynamic limit 



N~ I D{E)f{E)dE 
'o 



is calculated by 



la (l + a;2 -/3(gf- l)/i)i-5 
where — —1 + P{q — l)/i and 



—dx (100) 



K 



A-kV /2mkT 



\q-l 



3/2 



(101) 



This integral could be computed as a beta integral 
if 1 < g < 7/5 

T{u + v) 

(102) 



3(g- l)/i3 V 9 - 1 



where u — {5 ~ iq)/2{q ~ 1) and v — 1/(1 — g). 
Same calculation could be done for the generalised 
energy in the thermodynamic limit 



{2mEfl'^P{E)dE (103) 



after simple manipulation 

a-kV f2mkTV''^ _2«r(u)r(u) 



2/3((7- l)/i3 \q-l 



The average energy per particle 



= ^NkT 



(104) 



(105) 



It is worth noting that the generalized energy is in 
this equation related to the Tsallis constant (and 
not to the Boltzmann-Gibbs constant). 

7 Conclusions 

In this paper, we studied the probabilistic proper- 
ties within the framework of Tsallis thermodynam- 
ics. It has been shown that movariance and mari- 
ance are suitable functions for generalized covari- 
ance and variance in the non-extensive thermody- 
namic context. The invariance of thermodynamic 
relations is established if the Tsallis constant de- 
pend on the entropic parameter. However, in the 
canonical ensemble, as illustrated on the ideal gas, 
the equation of state as the same form as the ordi- 
nary equation in the non-extensive formalism with 
the Bolzmann-Gibbs constant. In the grand canon- 
ical ensemble, the dependance of the q-average in- 
ternal energy is related to the Tsallis constant. In 
both cases, we recover the usual equations in the 
extensive limit. 
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